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Abstract. We discuss necessary conditions for the existence of 
probability distribution on particle configurations in d-dimensions 
i.e. a point process, compatible with a specified density p and ra- 
dial distribution function <?(r). In d = 1 we give necessary and 
sufficient criteria on pg(r) for the existence of such a point pro- 
cess of renewal (Markov) type. We prove that these conditions 
are satisfied for the case g(r) = 0, r < D and g(r) = l,r > D, 
if and only if pD < e _1 : the maximum density obtainable from 
diluting a Poisson process. We then describe briefly necessary and 
sufficient conditions, valid in every dimension, for pg(r) to spec- 
ify a determinantal point process for which all ?i-particlc densities, 
jO n ( r ii ■■■) r n)j are given explicitly as determinants. We give several 
examples. 
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1. Introduction 



The microscopic structure of macroscopic systems, such as fluids, is 
best described by the joint n-particle densities p n (vi, ...r n ) where the 
ri, ...,r„ are position vectors in ^-dimensions p. The most impor- 
tant of these are the one particle density pi(ri) and the pair density 
P2( r i, r 2)- For spatially homogeneous systems, the only ones we shall 
consider here, From these one can also obtain an infinite set of condi- 
tions in the case where only p 1 and p^ are given . These conditions 
are very hard (or impossible) to check so the real question is whether 
one can get away with a smaller number of readily checkable conditions. 

A simple subset of such positivity conditions, emphasized by Percus 
[2] and by Stillinger, Torquato, et al. which follow directly from 

the definitions are, 



Conditions are obvious while ()1.2|) ensures that variances of one 
particle sum functions, ip — E^W) are non-negative, since it follows 
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from the definition of P2 that 



(1.3) 



<«'-«'>M^)Y 



dkl0(k)l 2 £(k) 



where 0(k) = J e lk r <f)(r)dr, and the averages are with respect to the 



probability distribution of the point process. Simple considerations, 
see jS], show that one should add to these (jl.lj) and (jl.2j) at least one 
further requirement: the variance V A of the number of particles N\ in 
a region A, which corresponds to </>(r) being the characteristic function 
of the region A, must be larger than 0(1 — 0), where 8 is the fractional 
part of the average number of particles in A. That is, if (N\) = k + 8, 
for k a non-negative integer, k — 0, 1, 2, then 



The bound (jl.4j) comes from the fact that N\ can only take non- 
negative integer values, see jB] and Appendix where a more general 
condition of type (jl.4|) is proven. 

A simple one dimensional example for which (jl.lj) and (jl.2j) are sat- 
isfied but fll.4j) is violated, is 



(1.4) 



V A = ((N A -(N A )) 2 ) > 0(1-0). 



(1.5) 




r > p 



r < p 



-i 



-i 
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A direct computation shows that S(k) > 0, with equality holding for 
k — 0, (see below), while the variance in an interval of length L, L > 
p~ x is equal to 1/5, which violates ([1.4)1 . whenever #(1 — #) > 1/5. It 
is possible that the condition (jl.4|) becomes less important in higher 
dimensions where the minimal variance will go to infinity as the domain 
grows. For spherical domains it will grow at least like the surface area 
P~2) [IB]- Note that by choosing the value of r beyond which g(r) = 1 
as slightly smaller than p~ l (jl.l|) and (jl.2|) would be satisfied but (jl.4|) 
would not for some L. 

It is however possible, especially for the type of g(r) considered in 
[I] [H] that (|l.lj) and (|1.2jl are enough to ensure realizability. These 
g(r) have a hard core exclusion, prohibiting the centers of two particles 
from coming closer than a certain distance D, i.e. g(r) = for r < D. 
In particlar, it was conjectured in [I] [E] that it is possible to find a 
point process with density p > and a g(r) of the form 

{0, r < D 
1, r > D 

as long as pv(D) < 2~ d , where v(D) is the volume of a ci-dimensional 
sphere of radius D/2. For pv(D) > 2~ d , S(k) will be negative for k = 0. 
(11.61) also satisfies condition ()1.4j) for pv(D) < 2~ d , although this was 
not explicitly imposed. There are also heuristic arguments, bolstered 
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by computer simulations [Zj and by considerations of the d — > oo limit 
[T4"j . for the realizability of (jl.6j) when pv(D) < 2~ d . 



since it yields a system for which the variance V\ grows only like the 
surface area of the boundary of A rather than the volume. Such systems 
are variously called superhomogeneous [TH], [I7j, or hyperuniform [S|. 
In d = 1, the variance in an interval of length L, Vl can actually be 
bounded uniformly in L as analyzed in ^3] , ^7] . Thus for the example 
ffOil with pD = 1/2, 



This is, by (jl.4|) . the minimum permissible variance when pL = k + 



Inspired by the work of Stillinger and Torquato we give here a proof of 
realizability of the model g(r) in f)l ,6|) for the case d = 1 and pv(D) = 
pD < e _1 . This is based on a particular construction of the point 
process as a dilution of a Poisson process with pD = XD exp[— XD]. It 
turns out that the new process is a Markov or renewal process [TH]- 
|19j . This permits us to describe all higher order correlation functions 
in terms of g(r). We do not know at present whether there exist non 
renewal point processes for some or all pD £ (e -1 , |]. We also do 



The case pv(D) = 2 d , for which S(0) = 0, is of particular interest, 



(1.7) 




2, k — 0, 1, 2, ... . 
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not know whether the explicitly constructed process for pD < e _1 is 
unique. In principle there can exist more than one process with the 
same p and g(r) but different higher order correlations; see sec. 5. 

We note that one dimensional renewal processes, described in sec. 
2, and determinantal processes for arbitrary dimension, described in 
sec. 4, are the only examples we know for which one can explicitly 
(and easily) construct higher order correlations from p\ and p^. In 
some cases these processes correspond to the distribution of particles 
in equilibrium systems. There is also a formula for the entropy of a 
renewal process in terms of g(r) [TT?] . 



2. Renewal Processes 

A translational invariant one dimensional particle system with den- 
sity p > 0, is described by a renewal process (RP) whenever the condi- 
tional probability density of finding a particle (or point) at a position 
q on the line, given the configuration of all particles to the left of q, 
say, q-\ < qo < q, depends only on x = q — q± [TH] . [T^j . Let us call 
that density P\(x). In other words, given that there is a particle at q, 
P\{x) is the probability density that the first particle to the right (left) 
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of q is at q + x (q — x). This corresponds, if we think of the points as 
events in time, to a Markov process. Clearly 



(2.1) / P 1 (x)dx = l, / xP(x)dx = p~ 1 
Jo Jo 

Calling P n (x) the probability density for finding the nth particle at 

a distance x to the right of the specified position of a given particle we 

have 

(2.2) P n (x)= f P n - 1 (x-y)P l (y)dy, n = 2,3,... 

Jo 

By the definition of pg(r) we have 

oo 

(2.3) pg (r) = J2Pn(r) 

n=l 

Taking the Laplace transform of ()2.3j) . using ([2.2)1 . then gives 

poo 00 

pg(s)=p e~ sr g{r)dr = Y}Pi{s)} n 

J ° n=l 



(2.4) =A(*)/[i-A00] 

Conversely, a given p and ^(r) will be realizable as a renewal point 
process if and only if 

(2.5) Q(s) = pg(s)/[l + pg(s)] 



s 
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is the Laplace transform of a probability density, P\{r) > 0, satisfying 
(I2.1jh We will show in the next section that for the one dimensional 
g(r) given in (|1.6|) this is true when and only when pD < e -1 . 

It is clear from the definition of a renewal process that the higher 
order correlation functions of such a system can be readily expressed in 
terms of p and g(r). More specifically given points X\ < x 2 < ... < x n 
on the line we have for n = 3, 4, 

(2.6) p n (xi, ...,x n )/p n _i(xi, ...,x n -i) = pg(x n -x n -i) 

since the left side is just the particle density at x n given that there are 
particles at x±, x n ^\. Thus 

(2.7) p 3 (x 1 ,x 2 ,x 3 ) = p 3 g(x 2 - x 1 )g{x 3 - x 2 ), x x < x 2 < x 3 , 



There is also a simple expression for s, the entropy per unit length 
of a renewal process [T^j, [Hj- It is given by the following formula, see 
Aizenman-Goldstein-Lebowitz 



Jo 

where Wo(^) = f°° Pi(y)dy is the probability that there is no particle 
between q and (q + r). 



etc. 



(2.8) 
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We can realize an RP as an equilibrium system of particles in d = 1 
in which only nearest neighbors interact: there are no interactions be- 
tween non-nearest neighbor particles, whatever the distances between 
them. For such a pair interactions u(r), P\(r) is given by P 

(2.9) P 1 {r) =Ce~ l3[pr+u{r)] , r > 0, 

where (3 is the reciprocal temperature, p = p{(5, p) is the pressure and 
^ = Uo° e~P' pr+u ( r "dr]~ 1 is a normalization constant. Conversely given 
P\(r) we can always define a f3u(r) and the corresponding /3p by in- 
verting $nty . 

A well known example of such an equilibrium system with only pair 
interactions is that of hard rods with diameter D, u(r) = oo for r < D, 
u(r) = for r > D. For this system P\{r) = 0, for r < D, and 

(2.10) Pi(r) = (3pe~ Mr - D \ for r > D 
with 

(2.11) (3p = p[l- P D}-\ 

Eq. (|2.7|) then gives the well known formula for the entropy density of 
this system [l] 

(2.12) s = _ p i og [ p /(i_ pjD )] +p 
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3. The Realizability of (jl.fij) as an RP 



By general theorems a necessary and sufficient condition for a func- 
tion of s, to be the Laplace transform of a non-negative density is that 
it be "completely monotone" for all s > [20]. That is, it is required 
that its derivatives alternate in sign for all s > 0. Thus for g to define 
a renewal process it is necessary and sufficient that Q(s) in ()2.5|) have 
the property that 



(3.1) (-l) fc Q (fc) (s) > 0, for all k = 0,1, 2, and all s > 0, 

where / (fc) (s) = d k f (s) / ds k . 
For the g(r) in ([l.bjl . 



It can be shown that (J3.1j) will be satisfied by ()3.3)1 if and only if 
pD < e _1 Here we provide a simple construction of this point 

process by starting with a Poisson process on the line, x G (— oo, oo), 
with density A and removing points which are too close ending up with 
a density p = \e~ XD and the step g{r) of (jl.6|) . 



(3.2) 




and the corresponding Q(s) in (|2.5j) is 



(3.3) 



Q(s)=pe- sD /[s + pe- sD ] 
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The procedure is as follows. Denote the points of the Poisson process, 
by — x 2 , — xi, x , Xi, x 2 , x 3 , with x, < x i+ \. Then if 
(x i+ i — x^ < D, Xi is removed; if (xj+i — Xj) > D, x, stays. Now 
the probability that (xi+i — Xj) is greater than D is e~ AD so the density 
of remaining points is 

(3.4) pD = XDe~ XD < e -1 . 

The last inequality follows from the fact that ye~ y has its maximum 
value e _1 at y — 1. Note that for pD < e~ l there are two different 
values of A which lead to the same RP with density p (see below). 

The new translation invariant process with density p clearly has 
g(r) = for r < D. To see that g(r) = 1 for r > D we note that, given 
a surviving point at position q, the density of other surviving points at 
q + r is, for r > D, just the density of points for the Poisson processes 
which have survived, i.e. Ae~ AD = p. 

It is clear from the above construction that the new process satisfies 
the conditions at the beginning of sec. 2 and so is an RP with P\(r) = 
Q(r), the inverse Laplace transform of Q(s) in ()3.3|) . To compute Q(r) 
we use units in which p — 1. Define 



(3.5) 



Q(r) = Q(y + nD) = w n (y), for nD < r < (n + l)D 
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and < y < D, n = 0, 1, 2, .... It is then easy to deduce from 
that 

rv 

(3.6) w n +i(y) = «Wi(0) - / w n (x)dx, n>0 
with 

(3.7) *uo(y) = 0, = 1, W2(y) = 1 — 3/, ... 
Furthermore 

(3.8) w n+1 (0) = w n {D) for n > 1 

i.e., Q{r) is continuous for r > D. 
Define now, 

oo 

(3.9) ^(A;y) = j> B «; B (y) 

n=l 

It follows then from (|3.fj|) that 

rv 

(3.10) i/>(\,y) = i/)(\;0)-\ ip{X;x)dx 

Jo 

and thus 

(3.11) V(A;2/) = ^(A;0)e-Y 
Putting y = D then yields 



(3.12) ^(A; -D) = ^(A; 0)e" AD = i[V>(A; 0) - A] 
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where the last equality follows from ()3.8|) and (J3.7|) . This gives 

(3.13) ^(A;0) = T -A_ J . 

The positivity of Q(r) is equivalent to the requirement that all the 
coefficients Cj in the expansion of 

oo 

(3.14) ^(A;0)=^C^ 

3=0 

are positive. This again leads to the requirement that D < e _1 , with 
the explicit formula (due to E. Speer) 

n 

(3.15) Wn (y) = J2l(n-k)D + y] k - 1 (-l) k - 1 /(k-l)\, D < e~ l . 

k=i 

4. Determinantal Point Process 

We review here briefly how one can obtain point processes from a 
g(r) satisfying certain inequalities in any dimension. The construction 
of such processes is a subject of great current interest in mathematics 
and we refer the reader to |22] f° r more information. We again restrict 
ourselves to homogeneous systems and choose units in which p = 1. 
Let B(r) be a complex function such that 



(4.1) 



B(r) =5*(-r), 5(0) = 1 
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< B(k) 



e- ikr B(r)dr < 1 



It can then be proven that conditions (|4.1|) and (|4.2|) are necessary and 
sufficient for the existence of a point process with n-particle densities 
given by the determinants [22] 



(4.3) 



p n [n, ...,r n ) 



1 B(r 12 ) ... B(r ln ) 



B(r 21 ) 1 ... B(r 2n ) 



B(r nl ) B(r n2 ) 



where = r-j — Yj. In particular we have 



(4.4) 



g(r) = 1 - \B(r) 



with g(0) = 0, g(r) < 1. 

Such a process is called a determinantal point process (DP). It follows 
that given a g(r), such that the Fourier transform of B(r) = [1— g(r)] 1 / 2 
satisfies (14. 2|) and ^(0) = 0, < g{r) < 1, we can construct a point 
process with explicit correlations (J4.3|) . This gives a large (uncountable) 
class of g(r) which have the realizability property. For all details, see 
and references there. 



We make two remarks: 
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1) To obtain a superhomogeneous system [H] [H! with ,§(0) = 0, for 
the determinant al point process specified by some B(r), it is necessary 
and sufficient that B(k) be a characteristic function of a set Q in k- 
space, i.e. -B(k) = 1 for k G Q, B(k) = otherwise. This is the 
case for the well known one dimensional system of particles on a circle 
with pair interaction 4>{rij) = — e 2 log | r^- 1 , at reciprocal temperature 
(3 = 2e~ 2 . For this system, with p = 1, the infinite volume limit of 
the radial distribution function is given by g(r) = 1 — (sin7rr/7rr) 2 and 
the variance Vl of the number of particles in an interval of length L 
grows like log L. This system is sometimes referred to as the Dyson gas: 
the p n describe the correlations of the eigenvalues of random Gaussian 
Hermitian Matrices |23|. 

2. To get translation invariant determinantal correlation functions 
as in ()4.3|) it is not necessary that B depend only on Y\2- It is only 
necessary that i?(ri,r2) = F(ri2)e^ <J?,( ' ri ' ) ~ < ^ r2 - ) l with i?(ri,r2) satisfying, 
as an operator, the analogue of (|4.1|) [22j, [23 • This is the case for a 
two dimensional one component plasma with 4>(rij) = — e 2 log |, (3 = 
2e~ 2 [13, [23 • For this system the variance in the number of particles 
in a disc of radius R grows like R [T3] . 
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5. Example and Discussion 

We illustrate here the construction of a DP in d dimensions from a 
given pi and p<i which is, in d — 1, also an RP. As in the example (jl.6|) 
this can be done for only a subset of the parameter for which (jl.lj) . 
fll.2j) and (jl.4j) are satisfied. On the other hand everything here can be 
computed explicitly in an elementary way. Using units in which p = 1, 
let 

(5.1) g( r ) = l-e- Xr , A>0. 

It is easily checked that this g satisfies (jl.lj) . (jl.2j) and (jl.4j) . whenever 
A > A d , Ai = 2, A 3 = (87T) 1 / 3 , .... It follows from (jOJ) that this g 
determines a DP with B = e _Ar//2 whenever A > 2\d- 

On the other hand, using (j2.5j) . we get for d = 1, the Laplace trans- 
form, 

(5.2) Q(s) = \/[s 2 + Xs + X] 

from which one readily finds, by factorizing the denominator in (j5.2j) 
and using criteria (j3.1j) . that (j5.1j) determines a RP if and only if A > 
2Ai = 4. For such values of A, it is then easily found that 
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p x (r) = A [A 2 - 4A] 




A + VA 2 - 4A 
2 



r 



exp[ 



A - VA 2 - 4A 
2 




A > 4 



(5.3a) 



Pi(r) = 4re 



-2r 



A = 4 



(5.3b) 



It is now easy to check that ()2.6j) and ()4.3|) give the same correlations 
p n (xi, ...,x n ). In particular for all A > 4, 



(5.4) p 3 (x 1 ,x 2 ,x 3 ) = (l-e- x ^-^)(l-e- x ^- x ^), Xl < x 2 < x 3 



Using ()5.3)1 in (|2.8j) one can also obtain the entropy of this system for 



The fact that the DP and RP constructions in the above example 
yield the same point process might suggest that an RP or DP deter- 
mines a unique point process. This may indeed be the case. We note, 
however, that uniqueness is not true in general, as can be seen from 
considerations of systems with non-reflection invariant correlations. 
Thus, while we always have, for translation invariant systems, that 
g(r) = g(—r), there is no such symmetry for the higher order p n . In par- 
ticular there is an explicit construction of a translation invariant point 
processes in d = 1 which, when run "backwards" will have the same p 
and g(r) as the original process but a p' 3 obtained from the original one 



A > 4. 
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by reflection, i.e. p' 3 (x 1 ,x 2 ,x 3 ) = p 3 (-xi, -x 2 , -x 3 ) ^ p 3 (x 1 ,x 2 ,x 3 ); 
see [12] for details. 

The sufficiency of conditions (|1.1|) . (|1.2|) and (jl.4j) (with the gen- 
eralizations given in the Appendix) remains open although it seems 
unlikely that any finite number of conditions would suffice for the gen- 
eral case 0,0, [HI- The construction of the step g{r) in (|1.6|) by the 
dilution of a Poisson process does not seem to work in d > 1. On the 
other hand we have not found any counterexample so far. 
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Appendix: Proof of (11.41) 

We give here an elementary proof of (jl.4|) and its generalizations, see 
also Let P{k) = Prob. of having k particles in A such that 

oo 

(k) = kP{k) = N + e, n = o, i, 2, o < e < 1. 

fc=0 
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Then the variance 

V A = J2[k-N-6] 2 P(k) 

k 

= J2( k - N) 2 P(k) - 26[{k) -N}+6 2 
= J2( k ~ NfP{k) - e 2 
> J2( k ~ N)P(k) - 6 2 = 6(1 - 6). 

The inequality follows from the fact that n 2 > n for n an integer. 
Equality occurs when P{k) = a5k,N + (1 — qO^.tv+i, with a determined 
by#. 

We note that the same argument works also for the variance of a 
linear combination of the number of particles in regions Aj. Let 
Y = Y2i=i m i^Ai, where rrii are integer coefficients, then ((Y — (Y)) 2 ) > 
6(1 — 0). In particular consider the difference between the number of 
particles in a region Ai and a region A 2 . Letting ((iV Al — N\ 2 )) = 
K + 9, K = 0, ±1, ±2, ... we again have 

V lt 2 = ([N Al -N M -K- 6} 2 ) > 6(1 - 6) 
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